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Abstract.  Local  feedback  stabilization  of  bifurcated  solution  branches  is  studied.  Two 
cases  are  considered;  that  in  which  the  nominal  system  undergoes  a  Hopf  bifurcation  as 
a  parameter  is  varied,  and  the  case  of  a  stationary  bifurcation  from  a  simple  zero  eigen¬ 
value.  For  each  case,  results  on  the  existence  of  a  stabilizing  feedback  are  given.  More¬ 
over,  simple  synthesis  techniques  for  the  stabilizing  controllers  are  discussed.  A  concept 
of  proximity  stabilization  is  introduced  as  an  alternative  to  stabilization  in  the  ordi¬ 
nary  sense  for  systems  that  are  not  locally  stabilizable.  A  result  is  stated  on  the  generi- 
city  of  proximity  stabilizability.  Motivation  for  further  research  in  several  areas  is 
given,  r  'x .  c  '  .  -■  L 


1.  Introduction 

A  standard  preliminary  step  in  the  analysis  and  design  of  control  systems  is  the 
linearization  of  the  model  dynamics  about  a  nominal  reference  trajectory.  The  analysis 
and/or  design  are  then  performed  for  the  resulting  approximate  linear  model.  The  suc¬ 
cess  of  this  technique  in  many  applications  can  be  attributed  to  the  result  of  Liapunov 
that,  if  the  linearized  system  is  locally  asymptotically  stable,  so  too  is  the  original  non¬ 
linear  model.  This  is  the  celebrated  “principle  of  linearized  stability."  This  principle 
bolds  for  finite  systems  of  ordinary  differential  equations,  as  well  as  for  some  infinite 
dimensional  problems.  Situations  in  which  the  local  stability  properties  of  the  nonlinear 
model  cannot  be  inferred  from  its  linearization  are  referred  to  as  critical  cases  in  stabil¬ 
ity  theory  [15]  . 
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Early  results  on  stability  in  critical  cases  were  obtained  by  Liapunov  [25]  ,  Malkin 
(26)  ,  Pliss  [33]  and  Krasovskii  [23]  .  This  topic  has  continued  to  be  an  active  area  of 
investigation.  Recently,  the  related  question  of  stabilization  in  critical  cases  has 
received  significant  attention.  The  papers  [12,6,35,4,1,2,16,7]  represent  a  broad 
spectrum  of  approaches  to  this  problem. 

Our  aim  in  this  lecture  is  to  summarize  some  recent  results  of  the  author  and  J.-H. 

Fu  [1,2]  on  the  closely  related  problem  of  local  bifurcation  control,  to  present  an  exten¬ 
sion  to  one  case  where  the  results  of  [1,  2]  do  not  apply,  and  to  indicate  some  directions 
for  further  work  in  this  area.  By  “local  bifurcation  control”  we  mean  the  local  stabiliza¬ 
tion  of  bifurcated  solutions.  Two  types  of  bifurcation  from  an  equilibrium  point  of  a 
one-parameter  family  of  ordinary  differential  equations  are  considered.  First,  the  stan¬ 
dard  Hopf  bifurcation  of  a  periodic  solution;  and  second,  stationary  (or  static)  bifurca¬ 
tion,  involving  only  equilibrium  points.  The  connection  between  local  stabilization  in 
critical  cases  and  local  bifurcation  control  becomes  transparent  given  some  basic  facts 
about  bifurcations  of  equilibria  of  differential  equations.  There  are  advantages  to 
analyzing  these  two  problems  simultaneously.  Often  results  on  stabilization  in  critical 
cases  can  be  directly  applied  also  to  problems  in  the  control  of  bifurcations. 
Correspondingly,  local  bifurcation  control  problems  provide  added  motivation  for  the 
study  of  stabilization  in  critical  cases. 

Our  approach  to  the  local  feedback  stabilization  problem  has  the  novel  feature 
that  it  facilitates  the  derivation  of  generally  valid  analytical  criteria  for  stabilizability, 
as  well  as  specific  stabilizing  feedback  controls.  This  is  possible  through  use  of  bifurca¬ 
tion  formulae  which  involve  only  Taylor  series  expansion  of  the  vector  field  and  eigen¬ 
vector  computations.  The  bifurcation  formulae  used  in  [1]  for  the  study  of  control  of 
Hopf  bifurcations  were  obtained  by  Howard  [19]  by  a  harmonic  balance  approach  using 
the  Fredholm  Alternative.  These  formulae  significantly  simplify  similar  formulae 
obtained  by  Hopf  [18]  .  In  our  study  of  stationary  bifurcations,  we  derive  bifurcation 
formulae  following  the  Projection  Method,  as  outlined  in  looss  and  Joseph  [21]  .  This  is 
essentially  the  same  technique  used  by  Howard  [19]  in  bis  study  of  Hopf  bifurcations. 

Previous  work  on  the  feedback  control  of  bifurcations  includes  that  of  Mehra  [28] 
and  Mehra,  Kessel  and  Carroll  [29]  .  See  also  the  account  in  Casti  [10]  .  These  studies 
tend  to  be  concerned  with  the  problem  of  globally  removing  bifurcations  by  state  feed- 
back.  The  general  results  apply  only  to  stationary  bifurcations,  since  they  are  obtained  /  > 

by  appealing  to  a  global  implicit  function  theorem.  This  is  in  contrast  to  the  local 
bifurcation  control  problems  considered  here,  where  one  seeks  only  to  modify  the  stabil-  6 
ity  properties  of  the  bifurcated  solutions. 

The  development  of  the  paper  is  as  follows.  The  problems  of  local  feedback  stabili-'  Faff 
zation  and  local  bifurcation  control  are  discussed  in  Section  2.  In  Section  3,  the  main  &i  df 
results  on  control  of  Hopf  bifurcations  are  recalled  from  [l]  without  proof.  In  Section  4  □ 

known  results  on  stationary  bifurcations  are  recalled,  and  bifurcation  formulae  are^*^  □ 

derived  for  this  case.  These  results  are  applied  to  the  stationary  bifurcation  control  “  °°  - 

problem  in  Section  5,  summarizing  the  main  results  of  [2]  .  A  notion  of  “proximity  sta¬ 
bilization”  is  introduced  in  Section  6,  and  shown  to  be  useful  in  problems  for  which 

stabilizing  control  law  exists  in  the  usual  sense.  Section  7  contains  two  worked.  -  -  - 
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examples,  one  illustrating  the  results  and  the  other  indicating  an  area  for  further  inves¬ 
tigation  related  to  the  proximity  stabilization  issue  of  Section  6. 

2.  Problem  Setting 

Consider  a  one-parameter  family  of  nonlinear  control  systems 

*=/<•(».«)  (1) 

where  z  E  R" ,  u  is  a  scalar  control,  /i  is  a  real-valued  parameter,  and  the  vector  field 
/  ^  is  sufficiently  smooth.  Suppose  that  for  ti  ^  0  Eq.  (1)  has  an  equilibrium  point  Xo(^) 
which  depends  smoothly  on  p.  In  the  sequel  the  system 

*=“/o(*.“).  (2) 

which  is  simply  (1)  with  p  =  0,  will  also  be  of  interest. 

This  paper  is  concerned  with  the  synthesis  of  feedback  controls  u=  u(i)  achiev¬ 
ing  certain  stability  properties  for  each  of  the  descriptions  (1)  and  (2). 

The  linearization  of  Eq.  (2)  at  z  =  0,  u  =  0  is  given  by 

z  =  Ax  ■+■  bu  (3) 

where  A  :=  and  6  :=  -^^O.O).  If  the  pair  {A  ,b)  is  controllable,’  then  a 

oz  au 

standard  linear  systems  result  |22]  asserts  the  existence  of  a  linear  feedback  u  =  -kx 
such  that  the  resulting  system  z  =  (^4  -  hit)z  is  asymptotically  stable.  Applying  this 
feedback  in  the  original  nonlinear  system  (2)  renders  the  origin  locally  asymptotically 
stable.  Moreover,  the  same  conclusion  applies  if  the  uncontrollable  modes  of  (3)  are 
asymptotically  stable.  In  contrast,  if  (A  ,6 )  has  an  unstable  uncontrollable  mode,  then 
the  origin  of  (2)  remains  unstable  regardless  of  the  applied  feedback.  These  considera¬ 
tions  imply  that  the  only  essentially  nonlinear  situation  encountered  in  local  feedback 
stabilization  occurs  when  some  uncontrollable  modes  of  (3)  are  pure  imaginary,  and  any 
other  uncontrollable  modes  are  asymptotically  stable.  However,  even  if  there  are  modes 
which  are  pure  imaginary  and  controllable,  it  may  be  important  to  study  the  existence 
of  nonlinear  stabilizing  feedback  control  laws  [34,  31,  32]  . 

It  is  natural  to  question  the  utility  of  results  giving  nonlinear  stabilizing  feedback 
controls  when  a  linear  stabilizing  control  law  exists.  There  are  indeed  several  reasons 
for  using  nonlinear  feedbacks,  two  of  which  are  noted  next.  First,  the  effect  of  a  linear 
feedback  control  designed  to  stabilize  the  linearized  version  of  the  critical  system  (2)  on 
the  one-parameter  family  of  systems  (1)  may  be  difficult  to  determine.  Indeed,  at  least 
for  small  feedback  gains,  one  can  expect  that  the  bifurcation  will  reappear  at  a  different 
value  of  the  parameter  p.  The  stability  of  this  new  bifurcation  is  not  easily  determined. 
Hence,  simply  using  a  linear  stabilizing  feedback  may  be  unacceptable  if  the  goal  is  to 
stabilize  a  bifurcation  and  not  merely  to  stabilize  an  equilibrium  point  for  a  fixed 

'Tkit  in**D*  tb»t  iht  itatt  •  of  Eq  (3)  can  b«  "atttrcd"  from  anj  iaitial  eoDdilioa  tg  to  aajr  prtacribtd  Urminal 
ceadilioD  witbia  a  pratpreiflfd  time,  b;  appropriat*  cbotct  of  tbr  eoatrol  a(i).  Siinilarl>,  oar  iprak;.  of  coatrollabilit)’ 
of  modt*  of  Eq  (3)  (ti|eo*alu«  dt  A),  dtptadiag  oa  wbtlbtr  or  aol  tbrar  modn  art  aStcltd  bp  alait  rttdback. 


parameter  value.  Second,  it  should  not  be  surprising  that  in  some  situations  a  linear 
feedback  which  locally  stabilizes  an  equilibrium  may  result  in  globally  unbounded 
behavior,  whereas  nonlinear  feedbacks  exist  which  stabilize  the  equilibrium  both  locally 
and  globally.  For  an  example,  see  Moon  and  Rand  [32]  .  Hence,  even  if  stabilization, 
rather  than  bifurcation  control,  is  the  issue  being  studied,  nonlinear  feedback  controls 
can  be  superior. 

In  the  light  of  the  foregoing  discussion,  it  is  appropriate  to  assume  that  the  matrix 
A  of  Eq.  (3)  possesses  at  least  one  eigenvalue  with  zero  real  part.  The  type  of  results 
obtained  will  depend  heavily  on  the  number  of  eigenvalues  of  A  which  are  assumed 
pure  imaginary,  their  multiplicity,  and  whether  they  are  zero  or  have  nonzero  imaginary 
parts.  The  results  of  this  paper  address  situations  in  which  either  of  the  following  two 
hypotheses  is  satisfied.  The  first  implies  that  Eq.  (1)  undergoes  a  Hopf  bifurcation  to 
periodic  solutions  at  =  0  when  u  =  0,  while  the  second  ensures  that  new  stationary 
solutions  of  (1)  bifurcate  from  Xo(/i)  ft=  0  when  u  =  0. 

(H)  Eq.  (1)  has  an  equilibrium  io(/')  ti  =  0.  Furthermore,  the  linearization  of 

(I)  near  Xo.  M  =  0  possesses  a  pair  of  simple,  complex  conjugate  eigenvalues  Xj(fj), 

^a(/^)  =  ^i(O)  =  >  0,  Re  X;  (0)  ^  0,  with  the  remaining  eigen¬ 

values  X3(0),  .  .  .  ,  X„(0)  in  the  open  left  half  complex  plane. 

(S)  Eq.  (1)  has  an  equilibrium  Xo(/i)  when  u  =  0.  Furthermore,  the  linearization  of 
(1)  near  x©,  a*  =  0  possesses  a  simple  eigenvalue  X,(/i)  with  Xi(0)  =  0,  X;  (0)  7^  0, 
with  the  remaining  eigenvalues  Xj^O),  .  .  .  ,  X„(0)  in  the  open  left  half  complex  plane. 

The  assumption  that  \j  (0)  ^  0  is  the  familiar  strict-crossing  (transversality)  condition 
introduced  by  Hopf  [18]  . 

Two  stabilization  problems  are  considered  in  the  sequel.  Both  are  studied 
separately  for  hypotheses  (H)  and  (S).  In  each  case,  one  of  these  problems  pertains  to 
Eq  (1)  and  the  other  to  Eq.  (2).  For  Eq.  (1),  the  goal  is  to  ensure  local  asymptotic  sta¬ 
bility  of  the  bifurcated  solutions.  This  will  be  referred  to  as  the  local  Hopf  bifurcation 
control  problem  or  the  local  stationary  bifurcation  control  problem,  depending  on  which 
hypothesis  is  in  force.  For  the  description  (2),  it  is  desired  to  solve  the  standard  local 
feedback  stabilization  problem  at  the  equilibrium  point  X(,(0).  Note  that  under  either 
hypothesis  (H)  or  (S),  Eq.  (2)  with  u  3  0  is  an  example  of  a  critical  nonlinear  system 
since  its  linearization  possesses  an  eigenvalue  with  zero  real  part. 

3.  Control  of  Hopf  bifurcations 

Under  hypothesis  (H),  the  Hopf  Bifurcation  Theorem  asserts  the  existence  of  a 
one-parameter  family  {  Pi  ,  0  <  <  <  <0  }  of  nonconstant  periodic  solutions  of  Eq  (1) 
emerging  from  x  *  0  at  =  0.  (This  assumes  ti  s  0,  of  course.)  Here  <  is  a  measure 
of  the  amplitude  of  the  periodic  solutions  and  cq  is  sufficiently  small.  The  periodic  solu¬ 
tions  p,{t )  have  period  near  and  occur  for  parameter  values  p  given  by  a  smooth 
function  |i(e).  Exactly  one  of  the  characteristic  exponents  of  p,  is  near  0,  and  it  is  given 
by  a  real,  smooth  and  even  function 
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^<)  =  -h  -h  •  •  •  (4) 

Moreover,  Pt{t)  is  orbitally  asymptotically  stable  with  asymptotic  phase  if  /9(e)  <  0  but 
is  unstable  if  /9(e)  >  0.  Denote  by  non  vanishing  coefficient  in  the  expan¬ 

sion  (4).  Checking  the  sign  of  P2K  sufficient  for  determining  stability.  Generically, 
K  =  1  so  that  locally  the  stability  of  the  bifurcated  periodic  solutions  p^  is  typically 
decided  by  the  sign  of  the  coefficient  $2- 

An  algorithm  for  the  computation  of  ^2  c&n  he  useful  in  the  solution  of  local  feed¬ 
back  stabilization  problems  under  hypothesis  (H).  In  [4]  the  evaluation  of  (a  scaled  ver¬ 
sion  of)  02  performed  using  a  formula  which  applies  to  two-dimensional  systems.  The 
original  n -dimensional  system  is  reduced  to  a  two-dimensional  system  by  appealing  to 
the  Center  Manifold  Theorem.  Use  is  then  made  of  the  fact  [8,  14,  27,  11]  that  the  sta¬ 
bility  properties  of  an  equilibrium  on  the  center  manifold  coincide  with  its  stability  in 
i?" .  In  fact,  the  value  of  0{t)  is  known  (27, 17,  3]  to  be  the  same  for  the  original  and 
the  reduced  systems.  The  approach  taken  in  this  work  (cf.  Abed  and  Fu  [l]  )  differs 
from  that  of  [4]  mainly  in  the  choice  of  algorithm  for  computing  02  The  implications 
for  the  type  of  results  one  obtains  are  nontrivial. 

Now  suppose  02  0.  Besides  locally  determining  the  stability  of  the  bifurcated 

periodic  solutions  P((f  )•  is  known  that  the  sign  of  the  coefficient  02  also  determines 
the  stability  of  the  equilibrium  xo{p)  at  criticality  {i.e.  at  /i  =  0).  This  fact  implies  that 
if  a  feedback  control  u  =  « (x )  can  be  found  such  that  02  <  0  for  the  Hopf  bifurcation 
occurring  in  the  controlled  system 

a  =/<«(»  .« (a ))  (5) 

then  the  local  feedback  stabilization  problem  described  in  Section  2  is  solved.  Simply 
use  the  feedback  u  ==  u(7)  in  Eq.  (1).  Indeed,  such  a  feedback  solves  both  the  local 
smooth  feedback  stabilization  problem  for  Eq.  (1)  and  the  local  Hopf  bifurcation  stabili¬ 
zation  problem  for  ony  parametrized  version  of  (1)  of  the  form  (5).  This  establishes  the 
connection  between  local  feedback  stabilization  and  Hopf  bifurcation  control. 

Rewrite  Eq.  (1)  in  the  series  form 

X  —  Lqx  u 'I  +  uLxX  -h  Qq(x,x) 

-h  Cffx  ,x,x)  -h  ■  ■  ■  (6) 

where  the  terms  not  written  explicitly  are  of  higher  order  in  x ,  u  and  p  than  those 
which  are.  Thus  Lq  and  Lj  are  square  matrices,  7  is  a  constant  vector,  Q^x,x)  is  a 
quadratic  form  generated  by  a  symmetric  bilinear  form  Qo{x  ,y)  giving  the  second  order 
(in  x)  terms  at  u  =  0,  /<  =  0,  and  Co(x,x,x)  is  a  cubic  form  generated  by  a  sym¬ 
metric  trilinear  form  C'o(x,y,x)  giving  the  third  order  (in  x)  terms  at  u  =  0,  p  =  0. 
(Note  that  Lo  is  simply  A  of  Eq.  (3),  and  7  corresponds  to  b  .) 

Denote  by  r  the  right  (column)  and  by  I  the  left  (row)  eigenvector  of  Lq  with 
eigenvalue  .  Normalize  by  setting  the  first  component  of  r  to  1  and  then  choose  / 
BO  that  ir  ss  1. 
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It  is  well  known  that  only  the  quadratic  and  cubic  terms  occurring  in  a  nonlinear 
system  undergoing  a  Hopf  bifurcation  influence  the  value  of  02-  Thus  only  the  linear, 
quadratic  and  cubic  terms  in  an  applied  feedback  v{x)  have  potential  for  influencing 
02-  To  simplify  the  analysis  and  to  emphasize  the  influence  of  nonlinear  terms  in  the 
feedback  control,  we  require  u  (x )  to  be  of  the  form 

«(x)  =  x''^,x  -#■  C,(x,x,x),  (7) 

where  is  a  real  symmetric  n  Xn  matrix,  and  is  a  cubic  form  generated  by  a 
scalar  valued  symmetric  trilinear  form. 

Theorems  1  and  2  below  give  sufi'icient  conditions  for  local  stabilizability  of  a  Hopf 
bifurcation  by  feedback  of  the  form  (7).  Both  are  positive  results.  Theorem  1  applies 
under  the  hypothesis  that  Z'}  ^  0,  while  in  Theorem  2  one  /'y  =  0  is  assumed.  By  the 
well  known  Popov-Belevitch-Hautus  (PBH)  eigenvector  test  for  controllability  [22]  of 
modes  of  linear  time-invariant  systems,  the  former  case  corresponds  to  the  critical 
modes  being  controllable  for  the  linearized  system  (3),  while  in  the  latter  case  these 
modes  are  uncontrollable.  The  theorems  were  proved  in  (ij  using  the  formula  for  02 
derived  by  Howard  (19]  .  Reference  [l]  also  contains  specific  formulae  for  stabilizing 
feedback  controls. 

Theorem  1.  Let  hypothesis  (H)  hold  and  assume  that  /-y^O.  That  is,  the  critical 
eigenvalues  are  controllable  for  the  linearized  system.  Then  there  is  a  smooth  feedback 
«(i)  with  u  (0)  =  0  which  solves  the  local  Hopf  bifurcation  control  problem  for  Eq.  (1) 
and  the  local  feedback  stabilization  problem  for  Eq.  (2).  Moreover,  this  can  be  accom¬ 
plished  with  only  third  order  terms  in  ti  (x ),  leaving  the  critical  eigenvalues  unaffected. 

Theorem  2.  Suppose  that  hypothesis  (H)  is  satisfied  and  that  Z'y  =  0.  Then  there  is  a 
smooth  feedback  u  (x )  with  u  (0)  =  0  which  solves  the  local  Hopf  bifurcation  control 
problem  for  Eq.  (1)  and  the  local  feedback  stabilization  problem  for  Eq.  (2)  provided  that 

0^Re{  -2/go(r.  jIoS) 

-f-  lQo{r  ,  j  {2iu},f  -  LoY^l) 

■h  ^l\2Lir  -h  I,r  )  }.  (8) 

4.  Stationary  bifurcations:  Analysis 

Under  the  stationary  bifurcation  hypothesis  (S),  it  is  well  known  [19]  that  Eq.  (1) 
exhibits  a  stationary  (or  static)  bifurcation  from  Xo  at  /i  =  0.  That  is,  new  stationary 
solutions  (i.e.  equilibrium  points)  bifurcate  from  Xq  at  =  0.  The  stability  characteris¬ 
tics  of  the  new  solutions  are  intimately  related  to  those  of  X(fp)  at  criticality,  i.e.  at 
—  0.  It  is  this  intrinsic  relationship  that  allows  the  joint  consideration  of  local  stabili¬ 
zation  for  Eq.  (2)  and  bifurcation  control  for  Eq.  (1). 

To  establish  this  relationship  and  motivate  the  derivations  to  follow,  consider  a 
general  one-parameter  family  of  nonlinear  ordinary  differential  equations 


(9) 

having  an  equilibrium  point  x^ii)  at  which  hypothesis  (S)  holds.  Then  near 
*o(®)>  0  in  (i ,  fi)  space  there  exists  a  locally  unique  curve  of  points  (  *(«), 

distinct  from  the  fi  axis  and  passing  through  (0,0),  such  that  for  all  sufficiently  small 
|e  ^  z(c)  is  an  equilibrium  point  of  (9)  when  fi  =/i(c).  Moreover,  the  parameter  t  may 
be  chosen  so  that  x{i)  and  /i(e)  are  smooth. 

Denote  the  series  expansions  of  /<(e),  x(e)  by 


(10) 

x(e)  =  Xje  -h  *2**  ■  ■  ■ 

(H) 

respectively.  Generically,  /ij  ^  0,  and  there  is  a  second  equilibrium  point  besides  io(p) 
for  all  small  ( /i|.  However,  if  /ij  =  0  and  /I2  >  0  (resp.  /12  <  0)>  there  are  two  new 
equilibrium  points,  one  for  positive  and  negative  values  of  e.  These  occur  only  for 
sufficiently  small  positive  (resp.  negative)  values  of  fi.  The  new  equilibrium  points  also 
have  an  eigenvalue  ^  which  vanishes  a.t  fi  =  0,  with  a  series  expansion 

0{t)  =  +  02^^  -h  •  •  •  (12) 

Moreover,  the  exchange  of  stability  formula  [19,  18] 

01  =  -  fti\j  (0)  (13) 

holds.  If  /ij  ==  0  and  fi2  0,  the  appropriate  exchange  of  stability  formula  is  [18] 

02=  ~2^2>^;(0).  (14) 

(Note:  Eqs.  (13)  and  (14)  may  be  derived  using  the  Factorization  Theorem  in  looss  and 
Joseph  [15,  pp.  90-91].)  Suppose  Xj  (0)  >  0.  Then  these  facts  imply  that  supercritical 
solution  branches  are  stable  while  subcritical  branches  are  unstable. 

The  following  result  follows  from  an  application  of  the  Center  Manifold  Theorem 
to  a  suspended  version  of  Eq.  (9)  at  Zo(0),  fi  =  0. 

Theorem  3.  Let  hypothesis  (S)  hold.  IJ  then  the  equilibrium  point  Zq{0)  is 

unstable  for  Eq.  (9).  //pj  =  0  and  P2  *o(0)  asymptotically  stable  if  /32  <  0 

but  is  unstable  if  02  >  0. 

Thus,  the  equilibrium  point  Xo(0)  will  be  assured  asymptotically  stable  if  one  can 
arrange  that  =  0  and  $2  explicit  formulae  can  be  derived  for  0i  and  02,  this 

provides  a  starting  point  for  the  construction  of  locally  stabilizing  feedback  controls  for 
the  critical  system  (2).  In  fact,  by  the  exchange  of  stability  formulae,  it  is  clear  that 
this  also  ensures  the  stability  of  the  bifurcated  stationary  solution,  by  ensuring  that  the 
bifurcation  is  a  supercritical  pitchfork  bifurcation.  This  is  a  desirable  outcome,  as  com¬ 
pared  to  the  transcritical  bifurcation  which  would  occur  if  0i  ^  0,  in  which  the  bifur¬ 
cated  equilibrium  point  is  stable  on  one  side  of  /i  =  0  and  unstable  on  the  other. 
Indeed,  under  hypothesis  (S),  a  supercritical  pitchfork  bifurcation  ensures  that,  even 
though  the  nominal  equilibrium  solution  Xo(p)  loses  stability  as  p  varies  through  0,  the 
new  equilibrium  solution  attracts  a  neighborhood  of  initial  conditions  about  Zo(0).  This 
can  also  be  shown  through  an  application  of  the  Center  Manifold  Theorem  and  the 


theory  of  normal  forms  [14]  for  details). 


Next,  bifurcation  formulae  for  Eq.  (9)  are  derived.  The  results  will  be  the  main  tool 
in  the  construction  of  stabilizing  feedbacks  for  Eqs.  (1)  and  (2)  in  Section  5.  The  Pro¬ 
jection  Method,  as  elaborated  in  [21]  ,  will  be  employed  in  the  derivation. 

By  assumption,  the  Jacobian  matrix  DgFfJ[xQ{0))  of  (9)  at  criticality  possesses  a 
simple  zero  eigenvalue  Xi(0).  Denote  by  r  (resp.  /)  the  right  column  (resp.  left  row) 
eigenvector  of  the  critical  Jacobian  matrix  corresponding  to  this  eigenvalue.  Using  the 
fact  that  0  is  a  simple  eigenvalue,  it  is  not  difficult  to  see  that  the  vectors  /  and  r  may 
be  chosen  to  have  only  real  elements.  To  be  more  specific,  set  the  first  component  of  r 
to  1  and  then  choose  /  so  that  Ir  =  1. 

Without  loss  of  generality,  assume  that  for  small  j  fi  \  the  known  equilibrium  point 
Zq{ii)  of  (9)  is  the  origin,  i.e.  XQ{fi)  =  0  for  small  |  fi  |.  This  can  always  be  achieved  by  a 
smooth  change  of  variables  z  —»  x  +  Xq{ii).  Rewrite  (9)  in  the  series  form 

x=L{n)x  +  Q^{i,x)  -f  C^{x,x,x)  -h  ••• 

—  Lqx  -h  uLix  +  -f-  •  ■  ■ 

->■  -f-  •  •  • 

•h  Co{x,x,x)  -f  •  •  •  (15) 

Here,  L{n),  £ j,  £2  ^re  n  Xn  matrices,  QJ[x,x),  Qo{z,z),  (?j(i,i)  are  vector  valued 
quadratic  forms  generated  by  symmetric  bilinear  forms  Q^{x,y),  Qo(^  ,y), 
respectively,  and  Cq{x  ,z  ,z  )  is  a  vector  valued  cubic  form  generated  by  a  symmetric  tri- 
linear  form  C{x  ,y,z).  The  terms  not  explicitly  written  in  (15)  are  of  higher  order  in  x 
and  fi  than  those  which  are. 

A  convenient  outcome  of  this  representation  is  the  formula 

X;(0)  =  U,,r.  (16) 

See  [19]  or  [21]  for  a  proof. 

If  I  is  any  real  (unknown)  solution  of  Fj{x)  =  0,  define  the  parameter  c  by 

(17) 

and  attempt  a  series  expansion  of  the  form 

( if') ) = .1.'*  ( ^". ) 

Substituting  the  expansion  (18)  in  the  equation  obtained  by  equating  the  right  side  of 
(15)  to  0,  and  equating  coefficients  of  like  powers  of  e  yields  the  following  relationships 

0  =  LoXi,  (19) 

0  =  Io»2  Qo(»i.  »i).  (20) 


.  4<.  41^ 


l.i  a.k  1 


Ik'  I4  .«t  kM  kU  . 


0  =  loZ3  ■>-  #‘1^1*2  -f-  /'2^1»1  lill'2^1 

+  2Qo{xi,Xi)  +  /iiQi(i,,i,)  +  Co(j,,T„x,).  (21) 

By  Eqs.  (17)  and  (18), 
c  =  lx  (0 

=  ilxi  +  1^1x2  +t^lx3  -h  ■•■  (22) 

Hence, 

Ixi  =  I,  and  Ix/i  =  0  for  k  >2.  (23) 

Eqs.  (19)  and  (23),  and  the  assumption  that  0  is  a  simple  eigenvalue  of  Lq,  now  imply 

*i  =  r  (24) 

Substituting  this  in  Eq.  (20)  gives  the  following  equation,  which  should  evidently  be 
solved  for  both  x^  and  pii 

^0^2=  -  Qd.rj)  (25) 

Recall  that  Lq  is  singular.  From  elementary  linear  algebra  (or  the  Fredholm  Alterna¬ 
tive),  this  equation  has  a  solution  x^  if  and  only  if  the  right  side  of  (25)  is  orthogonal  to 
all  left  eigenvectors  of  Lq  corresponding  to  the  zero  eigenvalue.  Since  zero  is  a  simple 
eigenvalue  of  Lq,  one  need  only  require  that 

/iitLjr  +  /Qo( '■.»■)  =  0,  (26) 

so  that  fii  is  determined  as 

where  Eq.  (16)  has  been  employed. 

Since  the  Fredholm  Alternative  conditions  are  now  satisfied,  Eqs.  (20)  and  (23)  for 
I2  have  a  solution.  This  solution  is  easily  verified  to  be  unique.  Equations  (20),  (23)  are 
conveniently  expressed  as  the  single  equation 


I ' J^“l  ®  I 


Since  (28)  has  a  unique  solution,  the  coefficient  matrix 


is  full  rank.  Hence,  R  is  a  nonsingular  square  matrix  and  X2  is  given  by 


=  (R’-Rr’«"  ( 0  ‘^'*'■'’1 


With  *2  now  available,  one  applies  the  Fredholm  Alternative  to  Eq  (21)  to  solve 


for  the  coefficient  Multiplying  both  sides  of  (21)  by  /  and  solving  for  pj.  one  obtains 

/‘2=  -  -r-i—  {  +  fiilLiT  -h  2lQo{r,X2) 

Aj  (Oj 

MilQi{r,r)  +  ^C’o(r  ,r  ,r  )  }. 

Using  the  exchange  of  stability  formula  (13)  and  Eq.  (27)  for  /ij,  the  coefficient 
is  found  to  be 

/?!  =  iQo{r  ,r  )  (32) 

If  /ij  =  0  (implying  also  =  0),  then  the  exchange  of  stability  formula  (14)  is  valid.  In 
that  case,  one  finds  that  given  by 

02  =  21  {2Qcir  ,X2)  -h  Co[r,r,r)),  (if  p,  =  0  only).  (33) 

The  formulae  (27),  (30)-(33)  will  be  employed  in  the  next  section  to  obtain 
sufficient  conditions  for  bifurcation  controllability  and  local  stabilizability  for  Eqs.  (1) 
and  (2),  respectively. 

5.  Stationary  bifurcations:  Control 

Motivated  by  the  bifurcation  formulae  derived  above,  and  by  the  results  of  Section 
3,  one  expands  the  vector  field  of  Eq.  (1)  as 

r  =  /  /t,u) 

=  LqX  -h  tiLiX  uL  ,i  -f-  u-j  -f  Qg(i,x) 

+  ft^Lsx  ^Q,(x,x)  -h  vQi(x,x) 

■h  Cf^x,x,x)  -h  •  •  •  (34) 

The  notation  here  is  similar  to  that  in  Eq.  (15).  As  in  Section  3,  a  feedback  control  con¬ 
sisting  of  quadratic  and  cubic  terms  is  assumed.  That  is,  u  =  u  (x )  is  taken  as 

u{x)  =  x^  QyX  C,(x,x,x),  (35) 

where  Q.  is  a  real  symmetric  n  Xn  matrix  and  C^{i,x,x)  is  a  cubic  form  generated 
by  a  scalar  valued  symmetric  trilinear  form.  Note  that,  as  in  Section  3,  u(i)  contains 
no  terms  linear  in  x .  This  ensures  that  the  left  and  right  eigenvectors  corresponding  to 
the  zero  eigenvalue,  and  the  value  of  ft  at  criticality,  will  be  unaffected  by  the  feedback 
control.  The  closed  loop  dynamics  with  a  feedback  of  the  form  (35)  become  (starred 
quantities  below  denote  values  after  feedback) 

x=ljx  -h  Ql{x,x)-h  Cl{x,x,x) 

-f  nLjX  -h  ft^Llx  -f  ^^;'(i,x)  -h  •  ••  (36) 

where  the  matrices  L*,  i  =0,1,2,  the  quadratic  forms  Qo(x,x),  Q*{x,i)  and  the 
cubic  form  Cl(x  ,x  ,x)  are 


I,-  =  li,  I  =  0,1,2,  (37a) 

+  Qo(*  .*)  (37b) 

<?;(»,*)=  (37c) 

and 

Co{z,x,x)=  C^(z,x,i)'r  +  Cq{x,x,z)  -h  (i *■  j )L  ,jr.  (37d) 

Symmetric  bilinear  and  trilinear  forms  Qoi^iV),  Cl{x  ,y  ,z)  generating  the  qua¬ 
dratic  and  cubic  forms  QHx  ,x)  and  Cl{x  ,x  ,z),  respectively,  are  now  chosen: 

Qo{^ ^y)  Quvh  -f-  Qoi^.y)-  (38) 


C*o{x,y,x)=  C^ix,y,2)'t  -h  C^x,y,z) 

—{{y^  (a^^C«y)^i^  -I-  (2^  Qui)L  ly  }■  (39) 

After  feedback,  the  coefficient  becomes,  using  Eq.  (32), 

=lQoir.r) 

—  i  {  Qo(r,r)  +  Q^rh) 

=  01  -f-  Q^r)l  -j,  (40) 

where  0i  denotes  the  value  of  0^  with  no  feedback,  i.e.  with  u  (j  )  s  0.  From  (40)  it  is 

clear  that  a  sufficient  condition  for  the  existence  of  a  feedback  « (jr )  driving  to  0  is 

hf^O.  (41) 

Indeed,  if  (41)  bolds,  then  any  feedback  control  of  the  form  (35)  with 

r^Qur=-p-  (42) 

ll 

results  in  =  0, 

Recall  from  the  Popov-Belevitch-Hautus  (PBH)  eigenvector  test  for  controllability 
of  modes  of  linear  time-invariant  systems  |22)  that  ^  0  is  equivalent  to  controllabil¬ 
ity  of  the  zero  eigenvalue  of  the  linearized  system  corresponding  to  Eq.  (2)  near  the  ori¬ 
gin.  Thus,  controllability  of  the  critical  zero  eigenvalue  for  the  linearized  system  is 
sufficient  for  the  existence  of  a  feedback  ensuring  0j  =  0.  Indeed,  if  ^  0,  this  condi¬ 
tion  is  also  necessary,  as  can  be  seen  from  Eq.  (40). 

As  outlined  in  Section  4,  the  next  step  after  ensuring  that  0*  —  0  is  to  arrange,  if 
possible,  that  0l  <  0.  By  Eq.  (33),  is  given  by 

0l  =21  {2Qlir,x;)  C;(r,r,r)).  (43) 

To  proceed,  it  is  necessary  to  evaluate  Zg,  according  to  the  formula  (30)  derived  in  Sec¬ 
tion  4.  Since  u(z)  contains  no  linear  terms,  the  matrix  R  occurring  in  Eq  (30)  is  the 


same  before  and  after  feedback,  as  is  clear  from  Eq.  (29).  The  fact  that  /i*  =  0  also 
simplifies  the  expression  for  x g  ■  One  has 

i;  =  I  j 

=  *2  -  {r^Q,r){R^R)-^R^  (2).  (44) 

From  (44),  it  is  clear  that  only  the  quadratic  terms  (i.e.,  Q^)  in  the  feedback  con¬ 
trol  influence  the  value  of  .  Using  (38)  and  (39),  one  evaluates  0£  : 

=21  {2^o(''.*f)  ^  2(r^Q^Tt)'r  -h  (r^C„r)Z,,r 
+  Co(r  ,r  ,r)  -h  C„(r,r,r)'r  } 

=  2/7  C,(r,r  ,r)  (45) 

where  ^  is  a  quantity  which  is  fixed  once  is  chosen,  and,  moreover,  does  not  depend 
on  the  trilinear  form  (z  ,y  ,z ). 

From  Eq.  (45)  it  follows  that  if  the  value  of  fig  may  be  assigned  arbi¬ 

trarily  by  appropriate  choice  of  the  scalar  valued  trilinear  form  C,  (x,y,z).  This  holds 
regardless  of  the  choice  made  for  the  symmetric  matrix  .  The  following  theorem  sum¬ 
marizes  these  results. 

Theorem  4.  Let  hypothesis  (S)  hold  and  assume  /^  ^^  0,  that  is,  the  critical  zero  eigen¬ 
value  is  controllable  for  the  linearized  version  of  Eg.  (2)  neor  the  origin.  Then  there  is  a 
amooth  feedback  control  u  =  u  (x )  with  u  (0)  =  0,  containing  only  quadratic  and  cubic 
terms  in  x ,  which  solves  the  local  stationary  bifurcation  control  problem  for  Eq.  (1)  and 
the  local  smooth  feedback  stabilization  problem  for  Eq.  (2).  Moreover,  the  quadratic 
terms  in  «  (z  )  can  be  used  to  ensure  that  =  0  for  the  controlled  system,  and  the  cubic 
terms  can  then  be  used  to  ensure  that  P2  <  0- 

Theorem  4  should  be  compared  with  Theorem  1,  which  contains  the  analogous 
results  for  local  Hopf  bifurcation  control.  The  assumption  ^  0  was  sufficient  for  ata- 
bilizability  in  that  setting  as  well.  However,  a  stabilizing  feedback  consisting  of  only 
cubic  terms  was  needed,  while  both  quadratic  and  cubic  terms  are  required  in  Theorem 
4.  This  is  due  to  the  need  for  a  two-stage  control  design  in  the  stationary  bifurcation 
control  case. 

Now  consider  the  case  /^  =  0,  i.e.  let  the  critical  (zero)  eigenvalue  be  uncontroll¬ 
able  for  the  linearized  system.  In  the  setting  of  Section  3,  under  the  analogous  assump¬ 
tion  it  was  found  that  generically  local  feedback  stabilization  of  the  nonlinear  system  is 
achievable.  However,  Eq.  (40)  reveals  that  in  the  present  setting  feedback  has  no  effect 
on  the  value  of  in  case  /^  =  0.  The  discussion  in  Section  4  therefore  implies  that 
the  local  feedback  stabilization  problem  for  Eq.  (2)  will  then  be  unsolvable,  unless 


perhaps  it  happens  that  =  0  in  the  absence  of  a  control  effort  (a  nongeneric  assump¬ 
tion).  Similarly,  the  local  stationary  bifurcation  control  problem  is  also  generically 
unsolvable  in  case  /“y  =  0. 

Theorem  5.  Let  hypothesis  (S)  hold  and  assume  l^  —  0,  that  is,  the  critical  zero  eigen¬ 
value  is  uncontrollable  for  the  linearized  version  o/ (2).  Then  if  ftk  0  for  Eg  (1)  with 
-u(z)  =  0,  both  the  local  stationary  bifurcation  control  problem  for  Eg  (1)  and  the  local 
feedback  stabilization  problem  for  Eg.  (2)  are  not  solvable  by  a  smooth  feedback  control 
with  vanishing  linear  part. 

Note  that  the  negative  conclusion  of  this  theorem  does  not  exclude  the  possibility 
that  “nearly  stabilizing”  feedback  controls  might  be  constructed  for  the  case  /-j  =  0 
In  the  next  section  we  consider  this  point,  and  introduce  precise  notions  of  “proximity 
stabilization”  and  “proximity  stabilizability”  as  an  alternative  to  local  stabilization  for 
this  case. 

6.  Proximity  stabilization 

Under  the  assumptions  that  a  simple  zero  eigenvalue  is  linearly  uncontrollable  and 
01  0,  Theorem  5  above  asserts  that  smooth  local  feedback  stabilization  is  unattain¬ 

able  within  the  class  of  purely  nonlinear  feedback  controls  This  fact  may  lead  one  to 
search  for  other  acceptable  forms  of  local  stability  besides  the  standard  notion  of 
asymptotic  stability  in  the  sense  of  Liapunov.  In  this  context,  recall  that  Brockett  |8’ 
has  obtained  easily  verifiable  necessary  conditions  for  stabilizability  by  any  smooth  feed¬ 
back  control.  Consider  the  following  definition. 

Definition  1.  The  origin  is  said  to  be  proximity  stabilizoble  for  Eq  (2)  if  for  any  <  >  0 
there  is  a  smooth  feedback  u  =  u  (r )  rendering  the  ball  centered  at  the  origin  of  radius 
e  in  a  locally  attracting  set. 

To  illustrate  the  nature  of  proximity  stabilization,  we  consider  a  simple  scalar  example. 
Suppose  X  satisfies 

1=1*-  «i.  (46) 

Then  01  0  and  /qf  =  0,  so  there  does  not  exist  a  smooth  feedback  « (x )  containing 

no  linear  terms  stabilizing  the  origin.  Indeed,  it  is  easily  checked  that  no  smooth  feed¬ 
back  can  render  the  origin  asymptotically  stable,  even  with  linear  terms  in  x  .  However, 
this  system  is  proximity  stabilizable  in  the  sense  of  Definition  1.  This  is  easily  seen  by 
noting  that  {  |x  I  <  t  }  is  attracting  for  (46)  if  u  (x )  =  ikx*  with  k  >e“*. 

The  following  result  states  that  under  the  assumptions  of  Theorem  5,  generically 
there  is  a  smooth  feedback  control  containing  only  cubic  terms  in  the  state  x  solving 
the  proximity  stabilization  problem.  Our  proof  is  not  difficult  but  in  the  interest  of 
brevity  it  will  be  presented  elsewhere. 

Theorem  8.  Under  the  hypotheses  of  Theorem  5,  the  proximity  stabilization  problem  for 
Eg.  (2)  is  solvable  with  a  cubic  feedback  control  under  generic  assumptions  on  the  func¬ 
tion  f  . 

7.  Examples  and  directions  for  further  research 
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The  first  example  considered  below  was  studied  previously  by  Su,  Meyer  and  Hunt 
[35]  .  They  derived  a  locally  stabilizing  feedback  control  law  using  feedback  linearization 
as  an  intermediate  step.  In  p>rinciple,  this  technique  allows  one  to  construct  an  infinite 
family  of  stabilizing  controllers  for  the  example.  This  example  is  one  of  a  nonlinear  criti¬ 
cal  system,  and  does  not  involve  a  bifurcation  parameter.  The  second  example  is  that 
of  a  parametrized  scalar  equation,  and  is  studied  mainly  to  motivate  consideration  of  an 
issue  of  proximity  stabilization  (in  a  sense  to  be  discussed  below)  of  bifurcations. 

Ejample  1.  Following  Su,  Meyer  and  Hunt  [35]  ,  consider  the  control  system 


|i;|  - 1 


sin  x> 
0  ‘ 


«(?)• 


(47) 


where  v  is  the  control.  Note  that  the  linearization  of  (47)  possesses  a  double  zero  eigen¬ 
value.  Thus  the  theory  of  Section  5  cannot  be  applied  directly.  To  circumvent  this 
problem,  we  define  a  new  control  u  by 


u  =  V 


-h  Q 

where  o  is  an  arbitrary  positive  number.  Then  (47)  becomes 
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(48) 


A  suitably  parametrized  version  of  (48)  would,  if  ii  =0,  undergo  a  stationary  bifurca¬ 
tion  from  a  simple  zero  eigenvalue.  Thus  the  results  of  Section  5  are  applicable  to  Eq. 
(48).  Note  that  no  quadratic  terms  appear  in  the  series  expansion  of  the  right  side  of 
(48)  for  u  =  0.  Hence  Q(,(7,y)s0.  By  (27),  this  implies  /ij  =  0,  and  hence  also 
^1  =  0  for  the  uncontrolled  system.  Theorem  4  then  implies  that,  if  /^r  ^  0,  only  cubic 
terms  would  be  needed  to  stabilize  the  origin  of  Eq.  (48).  Indeed,  it  is  easily  seen  that 
r  =  (1  0)^ ,  /  =  (1  q'^),  and  since  qr  =  (0  1)^ , 

/qr  =  q“*  ^  0.  (49) 

The  coefficient  ^2  uncontrolled  system  can  be  found  using  Eq.  (33).  We  have 

02  —  2/  C^T  ,r  ,r) 


=  2(1  Q->)(0) 


(50) 


Hence  stability  cannot  be  determined  based  only  on  knowledge  of  0i  and  02  However, 
using  Eq.  (45)  and  the  preceding  comments,  we  find  that  the  value  of  02  after  control  is 
applied  is 

0g  =  2  a~^  C,(r  ,r  ,r)  (51) 


and  we  are  free  to  choose  the  cubic  terms  C,  in  the  feedback  control.  (Recall  that  the 
quadratic  terms  have  been  set  to  0.)  It  is  easy  to  check  that  the  choice 

u  =  C,{x  ,x  ,x)  =  -  0.5  o  5  » 5. 


(52) 
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6  >  0,  renders  0l  <  0.  Hence,  for  the  original  system  (47),  we  have  the  stabilizing  con¬ 
trol  laws 

v  =  -Q-Xj  -  6  x]  (53) 

for  any  q,  6  >  0.  The  essential  distinction  between  this  control  law  and  that  obtained 
in  j35]  is  that  linear  feedback  in  only  one  variable  appears  in  our  control  law,  while 
linear  feedback  in  both  variables  forms  part  of  the  control  obtained  in  [35]  .  Computa¬ 
tions  involving  the  nonlinear  terms  are  essentia!  in  deriving  the  control  law  (53)  above. 

Example  £  (Motivational).  Consider  the  parametrized  scalar  equation 

X  =  /II  -  orx*  2o*x®  (54) 

where  p  is  the  bifurcation  parameter  and  o  is  a  design  parameter.  We  wish  to  view  (54) 
as  arising  from  a  specific  choice  of  the  form  of  a  feedback  controller  in  a  parametrized 
scalar  control  system.  The  term  2a®x*  has  been  retained  since  we  are  now  mainly 
interested  in  more  global  considerations.  It  is  easy  to  check  that  any  choice  of  a  >  0 
will  stabilize  the  (stationary)  bifurcation  occurring  at  /i  =  0,  x  =0. 

Suppose,  however,  that  we  are  also  interested  in  ensuring  the  persistence  of  some 
sort  of  local  asymptotic  stability  for  a  wide  range  of  values  of  the  parameter  p  >  0. 
Since  the  origin  is  unstable  for  /i  >  0,  we  will  be  satisfied  with  rendering  a  neighbor¬ 
hood  of  the  origin  attracting.  If  no  further  restrictions  are  placed  on  the  problem,  the 
solution  is  quite  simple:  given  an  interval  0  <  p  <  Pq  of  parameter  values  for  which 
attractivity  of  a  neighborhood  of  the  origin  is  desired,  there  is  an  a  sufficiently  large 
achieving  this  goal.  This  follows  by  a  simple  Liapunov  function  argument,  viewing  (54) 
as  a  perturbed  version  of  the  stable  equation 

a  =  -ox®.  (55) 

The  difficulty,  of  course,  is  that  the  size  of  the  neighborhood  which  is  attracting  may 

shrink  drastically  for  large  /i.  In  a  practical  sense,  therefore,  this  approach  may  prove 
unacceptable. 

Consider  the  following  alternative  formulation:  Given  a  “stability  region  tolerance” 
<0  >  0,  find  o  >  0  to  maximize  the  value  pi  such  that  the  set  [x  |  <  Cq  locally 
attracting  for  0  <  /i  <  p^.  The  following  analysis  addresses  this  issue.  Before  proceed¬ 
ing,  it  should  be  noted  that  the  important  goal  of  maximizing  the  size  of  the  region  of 
attraction  of  the  Cq  neighborhood  is  also  an  important  goal,  but  is  not  incorporated  into 
this  problem  formulation. 

Define  the  Liapunov  function  candidate 

V(i)  =  ,V2  (56) 

The  derivative  of  V  along  trajectories  of  (54)  is 

V^(x)  =  x®(/i  -  ox®  -f  2q®x*).  (57) 

Now  <  0  at  X  =  ±<0  if  and  only  if 
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11  -  acl  +  2  o’  <5  <  0, 

(68) 

i.e.,  if  and  only  if 

t*  <  Pin«(a)  :=  0€|  (  1  -  2  o  tj). 

(59) 

Now 

max,  >  0  ). 

(60) 

where 

o*  =  0.25  ef 

(61) 

is  the  best  value  of  a  that  could  be  obtained  with  the  Liapunov  function  V. 
tally,  fix  will  then  be  given  by 

Inciden- 

Pi  = 

(62) 

This  simple  example  is  given  as  motivation  for  studying  more  general  versions  of 
this  proximity  stabilization  issue  for  bifurcations.  The  basic  issues  noted  for  Example  2 
recur  in  the  analysis  of  the  general  case.  The  stability  of  critical  systems  depending  on 
parameters  has  received  much  attention  in  the  past.  Parameter  values  for  which  the  ori¬ 
gin  is  stable  for  the  critical  system  are  termed  safe  while  those  for  which  the  origin  is 
unstable  are  known  as  dangerous  [15,  5]  .  This  terminology  arose  from  the  consideration 
of  the  stability  properties  of  maneuvering  military  aircraft.  This  application  is  still  an 
active  area  of  investigation  [29,  30,  20}  .  The  author  is  continuing  these  efforts  by  inves¬ 
tigating  application  of  the  results  discussed  here  to  the  aircraft  high  angle  of  attack 
flight  control  problem. 
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Power  and  Energy  Systems,  Vol.  8,  No.  1,  pp.  37-43,  Jan.  1084. 

"Singularly  perturbed  Hopf  bifurcation,”  IEEE  Trans.  Circuits 
and  Systems,  Vol.  CAS-32,  No.  12,  pp.  1270-1280,  December  1085. 

“Multiparameter  singular  perturbation  problems:  Iterative  expansions 
and  asymptotic  stability,”  Systems  and  Control  Letters,  Vol.  5, 

No.  4,  pp.  270-282,  February  1085. 

“Decomposition  and  stability  of  multiparameter  singular  perturbation 
problems,”  IEEE  Trans.  Automatic  Control,  to  appear,  1086. 

(with  J.-H.  Fu),  “Local  feedback  stabilization  and  bifurcation  control, 

I.  Hopf  bifurcation,”  Systems  and  Control  Letters,  to  appear,  1086. 

“A  simple  proof  of  stability  on  the  center  manifold  for  Hopf 
bifurcation,”  SIAM  Review,  to  appear,  1086. 

(with  A.L.  Tits),  “On  the  stability  of  multiple  time-scale  systems,” 

Int.  J.  Control,  to  appear,  1086. 


STUART  S.  ANTMAN 


Professor  of  Mathematics 


Stuart  Antman  was  born  in  Brooklyn,  New  York  on  June  2,  1039.  He  received  his  B.S.  from 
the  Rensselaer  Polytechnic  Institute  in  1861.  He  received  both  the  M.S.  and  Ph.D.  degrees  from 
the  University  of  Minnesota  in  1863  and  1865  respectively. 

Before  coming  to  the  University  of  Maryland  in  1872,  Dr.  Antman  was  Assistant  and  Associ¬ 
ate  Professor  at  the  New  York  University  from  1867-1872.  He  has  held  numerous  visiting  profes¬ 
sor  positions  internationally.  Dr.  Antman  is  a  member  of  AMS,  MAA,  Pi  Mu  epsilon,  SIAM, 
Society  of  Natural  Philosophy  and  the  American  Academy  of  Mechanics.  Additionally,  he  partici¬ 
pates  in  many  activities  of  service  to  the  Academic  Community. 

RESEARCH  INTERESTS 

Professor  Antman *s  main  research  Interests  are  in  the  areas  of  Nonlinear  Elastostatisics  and  Non¬ 
linear  Bifurcation  Theory. 


SELECTED  RECENT  PAPERS 

(1)  "Boundary  Layer  Effects  for  Inflated  Elastic  Rings  and  Shells,"  with  M.C.  Calderer,  in 
preparation. 

(2)  “Beck’s  Problem  for  Nonlinearly  Viscoelastic  Rods,"  with  M.S.  Chen,  in  preparation. 

(3)  "Nonlinear  Problems  of  Elasticity,  Lectures  Notes,"  in  preparation,  400  pp. 

(4)  "The  Paradoxical  Asymptotic  Status  of  Massless  Springs,"  in  preparation. 

(5)  “  The  Singular  Global  Buckling  Problem  for  Anisotropic  Plates,"  with  P.  Negron,  in 
preparation. 

(6)  "  Travelling  Waves  in  Nonlinearly  Visocoelastic  Media  and  Shock  Structure  in  Elastic 
Media,"  with  R.  Malek-Madani,  in  preparation. 

(7)  "  Steady-State  Problems  of  nonlinear  Electro-Magneto- Thermo-Elasticity, 

(8)  "  Qualitative  Properties  of  Large  Buckled  States  of  Spherical  Shells,"  with  K.G.  Shih, 
University  of  Maryland  TR85-30,  61  pp.  To  appear  in  Archive  for  Rational  Mechanics  and 
Analysis. 

(8)  "Waves  in  Nonlinearly  Viscoelastic  Media.in  Nonstrictly  Hyperbolic  Laws."  with  R.  Malek- 
Madani,  Contemporary  Mathematics  Series,  AMS,  edited  by  B.L.  Keyflti  &.  H.C.  Kramer, 
to  appear. 

(10)  "The  Drawing  and  Whirling  of  Strings:  Singular  Global  Multiparameter  Bifurcation  Prob¬ 
lems,"  with  M.  Reeken,  University  of  Maryland  Technical  Report  1885,  58  pp.  SIAM  Jour¬ 
nal  of  Mathematical  Analysis,  to  appear. 

(11)  "Bifurcation  Problems  Associated  with  Nonlinear  Wave  Propagation,"  in  Proceedings  of 
the  Eighth  Dundee  Conference  on  Differential  Equations,  edited  by  B.D.  Sleeman  and  R.J. 
Jarvis,  Springer  Lecture  Notes  in  Mathematics,  to  appear. 

(12)  "Asymptotic  Shapes  of  Inflated  Spheriodal  Elastic  Shells,"  with  M.C.  Calderer,  Mathemati¬ 
cal  Proceedings  of  the  Cambridge  Philosophical  Society.  87  (1885),  541-548. 

(13)  "Asymptotic  Shapes  of  Inflated  Noncircular  Elastic  Rings”  with  M.C.  Calderer.  Mathemati¬ 
cal  Proceedings  of  the  Cambridge  Philosophical  Society.  87  (1085),  357-370. 
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JOHN  S.  BARAS 

Professor  of  Electricsl  Engineering 
Director,  Systems  Research  Center 


John  S.  Baras  was  born  in  Piraeus,  Greece,  on  March  13,  1048.  He  received  a  B.S.  in  Electr¬ 
ical  Engineering  with  highest  distinction  from  the  National  Technical  University  of  Athens, 
Greece,  In  1970.  He  received  the  M.S.  and  Ph.D.  degrees  in  Applied  Mathematics  from  Harvard 
University,  Cambridge.  Massachusetts  in  1071  and  1973  respectively. 

From  1073  to  1077,  1077  to  1082,  he  was  an  Assistant,  Associate  Professor  and  since  1082  a 
Professor  of  Electrical  Engineering  and  member  of  the  Applied  Mathematics  Faculty,  University 
of  Maryland.  Since  1085  he  has  been  Director  of  the  Systems  Research  Center.  He  has  held  visit¬ 
ing  research  scholar  positions  with  Stanford,  MIT,  Harvard  and  the  University  of  California. 
Berkeley.  Among  his  awards  are:  a  1078  Naval  Research  Laboratory  Research  Publication 
Award,  the  1080  Outstanding  Paper  Award  of  the  IEEE  Control  Systems  Society  and  a  1083  Alan 
Berman  Research  Publication  Award  from  NRL.  He  is  a  member  of  the  IEEE  Engineering  R  £  D 
Committee.  Professor  Baras  is  a  Fellow  of  IEEE. 

RESEARCH  INTERESTS 


Professor  Baras*  main  research  interests  are  in  stochastic  systems,  real-time  parallel  archi¬ 
tectures  for  nonlinear  signal  processing,  expert  systems  for  control  systems  synthesis,  distributed 
systems,  queueing  systems  and  scheduling. 


SELECTED  RECENT  PAPERS 

(1)  **  Complex  Variable  Methods  in  Modelling  and  Control  of  Linear  Distributed  Parameter  Sys¬ 
tems."  invited  paper,  to  appear  in  Frontiers  In  Applied  Mathematics. 

(2)  "  Decomposition  and  Decentralized  Control  System  Design;  A  Review  of  Frequency  Domain 
Methods,"  with  W.H.  Bennett,  Proc.  24th  IEEE  CDC,  1085,  pp.  1828-1835. 

(3)  "  A  New  Geometric  Stability  Margin  for  Feedback  Systems,”  with  W.H.  Bennett,  Proc. 
24th  IEEE  CDC,  1085,  pp.S10-S26. 

(4)  "  An  Impulse  Control  Problem  for  a  Stochastic  PDE  Arising  in  Nonlinear  Filtering."  with 
A.  Bensoussan,  Proc.  24th  IEEE  CDC,  1085,  pp. 1185-1187. 

(5)  "  Optimal  Sensor  Scheduling  In  Nonlinear  Filtering  of  Diffusion  Processes,"  with  a.  Bensous¬ 
san,  submitted  to  SIAM  J.  on  Control  and  Optimization. 

(6)  "  Expert  Systems  and  VLSI  Architectures  In  Real-Time  Nongausslan  Signal  Processing," 
with  A.  LaVigna,  Proc.  of  3rd  Bad  Honnef  Conference,  1085,  to  appear. 

(7)  "  Optimization  Based,  Direct,  Frequency  Domain  Design  of  Controllers  for  Linear  Distri¬ 
buted  Systems,"  1085  IFAC  Workshop  on  Control  Applications  of  Nonlinear  Program¬ 
ming  and  Optimization.  Pergammon  Press,  1086. 

(8)  "  Two  Competing  Queues  with  Linear  Costs  and  Exponential  Service  Requirements;  The 
/<c-rule  is  Often  Optimal,"  with  A.J.Dorsey  and  A.  Makowski,  Advances  In  Applied  Proba¬ 
bility,  Vol.  16,  No.  1,  (1084)  pg.8. 

(0)  "  K  Competing  Queues  with  Geometric  Service  Requirements  and  Linear  Costs;  Th  pc-ruie 

is  Always  Optimal,"  with  D.J.  Ma  and  A.  Makowski.  Systems  and  Control  Letters,  Vol.  6, 
No.  3,  (1085),  pg.  173. 


4 


CARLOS  A.  BERENSTEIN 


Professor  of  Mathematics 

Carlos  A.  Berenstein  was  born  In  Buenos  Aires,  Argentina  on  August  7,  1044.  He  received  a 
degree  in  mathematics  (Llcenciado)  at  the  Universldad  Nacional  de  Buenos  Aires  in  1066,  an  M.S. 
(1060)  and  the  Ph.D.  (1070)  from  the  Courant  Institute  of  Mathematical  Sciences  (New  York 
University).  In  1070-73  he  was  Benjamin  Pierce  Assistant  Professor  of  Mathematics  at  Harvard 
University,  since  that  time  he  has  been  at  the  University  of  Maryland  where  he  became  Associate 
Professor  In  1076  and  Full  Professor  In  1080.  He  has  been  a  Visiting  Professor  at  several  institu¬ 
tions:  Brandels  University,  University  of  Paris,  University  of  Bordeaux,  Scuola  Normale  Superiore 
(Pisa),  etc. 

Dr.  Berenstein  Is  a  member  of  the  American  Mathematical  Society.  He  is  currently  the 
organizer  of  a  Special  Year  In  Complex  Analysis  at  the  University  of  Maryland,  funded  by  the 
Nation2d  Research  Foundation  and  the  Argonne  University  Association. 


RESEARCH  INTERESTS 

Dr.  Berenstein’s  interests  He  in  Complex  Analysis,  Harmonic  Analysis  and  their  applications 
to  Engineering  These  applications  are  based  on  recent  developments  in  the  theory  of  functions  of 
several  complex  variables.  They  lie  in  Image  processing,  signal  processing,  frequency  domain 
methods  for  the  control  of  distributed  parameter  systems,  etc.  Some  of  his  students  are  currently 
working  on  the  Implementation  of  algorithms  based  on  his  work  using  MACS^'MA,  DELIGHT 
and  MAR\XIN. 


SELECTED  RECENT  PAPERS 

(1)  “  An  Inverse  Neumann  problem,"  with  P.  Yang,  to  appear  (1086).  33p. 

(2)  "  A  local  version  of  the  two-circles  theorem,"  with  R.  Gay,  to  appear  in  the  Israel  J.  Math, 
37p.  "  1-Inverses  for  polynomla  matrices  of  non-constant  rank,"  with  D.  Slruppa,  Systems 
and  Control  Letters,  Vol.  6  (1086),  pp. 309-314. 

(3)  "  Solutions  of  convolution  equations  in  convex  sets,"  with  D.  Struppa,  to  appear  in  Amer.  J. 
Math.,  37p. 

(4)  "  Spectra]  synthesis  on  symmetric  spaces,”  to  appear  in  Contemporary  Mathematics. 

(5)  "  Sur  la  synthese  spectrale  dans  les  espaces  symmetrlques,”  with  R.  Gay,  to  appear  in  J. 
Anal.  Math. 

(6)  "  Deconvolution  methods  for  multi-sensors,"  with  P.S.  Krishnaprasad  and  B.A.  Taylor, 
(1084)  67p. 

(7)  “  Ideal  generated  by  exponential-polynomials.”  with  A.  Yger.  to  appear  In  Advances  in 
Mathematics,  132  p. 

(8)  "  On  the  equilibria  of  rigid  spacecraft  with  rotors,”  with  P.S.  Krishnaprasad,  Systems  and 
Control  Letters.  Vol.  4  (1084),  pp.157-163. 

(0)  "  Le  probleme  de  le  deconvolution."  with  A.  Yger,  J.  Funct.  Anal.,  Vol.  54  (1083).  pp.ll3- 

160. 
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JOHN  HENRY  MADDOCKS 


Assistant  Professor  of  Mathematics 

John  Maddocks  was  born  in  the  United  Kingdom  in  1058.  He  received  his  B.S.  in 
Mathematics  with  honors  in  1078  from  the  University  of  Glasgow.  In  1082  he  received  his  Ph.D. 
from  Balliol  College.  University  of  Oxford. 

Before  becoming  an  Assistant  Professor  of  Mathematics  at  the  University  of  Maryland  in 
1084,  Dr.  Maddocks  was  an  Assistant  Professor  at  Stanford  University  (1081-82),  a  Junior 
Research  Fellow  at  the  Queen's  College,  University  of  Oxford  (1083-84)  and  a  Senior  Fellow  at 
the  Institute  for  Mathematics  and  Its  Application,  University  of  Minnesota  (1084-85).  He  is  a 
member  of  AMS,  SIAM  and  Society  for  Natural  Philosophy.  Additionally,  Dr.  Simo  participates 
as  a  reviewer  for  N.S.F.,  International  J.  of  Solids  and  Structures  and  J.  of  Applied  Math,  and 
Optimization. 

RESEARCH  INTERESTS 

Dr.  Simo's  research  interests  are  in  the  area  of  Elastostatlstlcs  and  Nonlinear  Bifurcation  Theory. 

SELECTED  RECENT  PAPERS 

(1)  "Stability  of  nonlinearly  elastic  rods,”  Arch.  Rat.  Mech.  Anal.,  85,  1084,  pp.  311-354. 

(2)  "Dynamical  theory  of  the  elastica,"  with  R.E.  Caflisch,  Proc.  Roy.  Soc.  Edin.  A.,  Vol.  00, 
pp.  1-23,  1084. 

(3)  "Restricted  quadratic  forms  and  their  application  to  bifurcation  and  stability  in  constrained 
variational  principles,"  SIAM  J.  Math.  Anal.,  Vol.  16  1085,  pp.  47-68. 

(4)  "A  model  for  twinning,"  with  G.P.  Parry,  J.  Elasticity,  in  press.  IMA  Preprint  Series  No. 
125,  1085. 

(5)  "Stability  and  folds,”  IMA  Preprint  Series  No.  174,  1085.  Submitted  to  Arch.  Rat,  Mech. 
Anal. 

(6)  "Mechanics  of  ropes,”  with  J.B.  Keller,  IMA  Preprint  Series  No.  183,  1085.  Submitted  to 
SIAM  J.  Appl.  Math. 

(7)  "Stability  of  the  elastica,"  Proc.  of  Nato/London  Math.  Soc.  Advanced  Study  Institute. 
Ed.  J.M.  Ball.  Reidel,  1083. 


JERROLD  E.  MARSDEN 

Professor  of  Mathematics 
University  of  California,  Berkeley 


Jerrold  Marsden  was  born  in  British  Columbia.  Canada  August  17,  1642.  He  received  his 
B.S.  from  the  University  of  Toronto  in  1865  and  his  Ph.D  from  Princeton  University  in  1667. 

Dr.  Marsden  has  served  as  an  Instructor  at  Princeton  University  (1867-68);  Lecturer  in 
Mathematics  at  Berkeley  (1668-66);  Assistant  Professor,  University  of  Toronto  (1870-71);  Assis¬ 
tant  Professor,  University  of  California,  Berkeley  (1871-72);  and  Associate  Professor,  University  of 
California,  Berkeley  (1872-76).  Since  1676  he  has  been  a  Professor  of  Mathematics  at  Berkeley. 
Additionally,  Dr.  Marsden  has  held  many  visiting  research,  professor,  fellow  and  scholar  titles 
nationally  and  internationally.  Dr.  Marsden  has  presented  numerous  invited  talks  and  conference 
lectures  and  received  a  number  of  awards  for  essays  and  research.  He  serves  the  scientific  com¬ 
munity  as  an  editor  or  reviewer  for  many  scholarly  publications.  Dr.  Marsden  has  written  14 
books  and  over  100  articles  on  various  topics  in  mathematics. 

RESEARCH  INTERESTS 

Dr.  Marsden's  research  Interests  are  in  the  areas  of  Nonlinear  Dynamics.  Hamiltonian  Systems. 
Chaotic  Dynamics,  Dynamical  Systems  and  Stability  among  others. 

SELECTED  RECENT  PAPERS 

(1)  "Constraints  and  mementum  maps  for  relativistic  classical  fields,”  with  Gotay,  Isenberg, 
Montgomery,  Sniatycki  and  Yasskin,  in  preparation. 

(2)  “Bifurcation  to  quasl-peridlc  Tori  in  the  interaction  of  steady  stale  and  Hopf  bifurcations," 
with  J.  Scheurie,  accepted,  SIAM  J.  Math.  Ann. 

(3)  "On  the  rotated  stress  tensor  and  the  material  version  of  the  Doyle-Ericksen  formula,”  with 
J.C.  Simo.  accepted.  Arch.  Rat.  Mech.  Ann. 

(4)  "The  York  Map  is  a  canonical  transformation,"  with  J.  Isenberg,  accepted.  J.  of  Geometry 
and  Physics. 

(5)  "Temporal  and  spatial  chaos  in  a  van  der  Waals  fluid  due  to  periodic  thermal  fluctuations,” 
with  M.  Slemrod,  accepted,  Adv.  in  Appl.  Math. 

(6)  "Averaging,  horseshoes,  and  exponentially  small  Melnikov  functions,"  with  P.  Holmes  and 
J.  Scheurie,  in  preparation. 

(7)  "Quasiconvexity,  second  variations  and  the  energy  criterion  in  nonlinear  elasticity,"  with  J. 
Ball,  accepted.  Arch.  Rat.  Mech.  Ann. 

(8)  "Nonlinear  stability  of  fluid  and  plasma  systems,"  with  D.  Holm,  T.  Ratiu  and  A.  Wein¬ 
stein,  accepted.  Physics  Reports. 

(8)  "The  Hamiltonian  structure  of  general  relativistic  perfect  fluids,”  with  Bao  and  Walton, 
preprint. 

(10)  "Generalized  Poisson  brackets  and  nonlinear  Liapunov  stability  •  application  to  reduced 
MHD"  with  Hazelline.  Holm  and  Morrison,  accepted,  Proc.  Plasma  Phys.  Conference, 
Lausanne,  June  1684. 

(11)  "Chaos  In  dynamical  systems  by  the  Polncare-Melnlkov-Arnold  method"  Proc.  ARO 
Workshop,  March  1684. 


JUAN  C.  SIMO 


Assistant  Professor,  Mechanical  Engineering 
Stanford  University,  California 


Juan  Simo  was  bom  in  Spain.  He  received  his  degree  of  E.T.S;  Ing.  Caminos  Canales  y 
Puertos  in  1976  from  the  Universidad  Politecnica,  Madrid.  Spain.  In  197C,  Dr.  Simo  received  his 
M.S.  degree  in  Metodos  Cuantitativos  from  Escuela  de  Organizacion  Industrial.  Madrid.  Spain. 
The  M.S.  and  Ph.D.  were  received  from  the  University  of  California  at  BerKely  In  Civil  Engineer¬ 
ing  in  1980  and  1982  respectively. 

Dr.  Simo  has  held  positions  as  Research  Assistant.  Lecturer.  Post-Doctoral  Fellow  and  Assis¬ 
tant  Research  Engineer  at  the  University  of  Calfironla  at  Berkeley  (1980-108S).  Most  recently  he 
has  been  at  Standford  University  as  an  Instructor  (1984).  Lecturer  (1986)  and  as  an  Assistant  Pro¬ 
fessor  (1985-present).  He  has  been  invited  to  present  his  work  at  a  number  of  conferences,  sympo¬ 
sium  and  lectures  in. the  United  States  as  well  as  internationally. 

RESEARCH  INTERESTS 

Dr.  Simo's  research  interests  are  in  the  areas  of  Nonlinear  Theory  and  Elastostatics. 


SELECTED  RECENT  PAPERS 

(1)  "Remarks  on  the  Patch  Test  and  Stummel’s  Counterexample,"  with  with  R.L.  Taylor, 
International  Journal  of  Numerical  Methods  In  Engineering,  to  be  submitted. 

(2)  "On  a  fully  Three-dimensional  finite  strain  viscoelastic  damage  model:  Formulation  and 
computational  aspects,"  Comp.  Meth.  Appl.  Mech.  Engng.,  submitted  for  publication. 

(3)  "On  the  Dynamics  of  Flexible  Beams  Subject  to  Large  Overall  Motions-  The  Plane  Case,” 
with  L.  Vu-Quoc,  Electronics  Research  Laboratory  Memorandum  No.  UCB/ERL  M85/63, 
University  of  California,  Berkeley,  submitted  for  publication  in  J.  Applied  Mechanics. 

(4)  "A  consistent  return  mapping  algorithm  for  plane  strain  elastoplasticity,”  with  R.L.  Taylor, 
Int.  J.  Num.  Methd.  Engng.,  (in  press). 

(5)  "A  3-Dimensional  Finite  Strain  Rod  Model:  Geometric  and  Computational  Aspects,”  with 
L.V.  Quoc,  Comp.  Mth.  Appl.  Mech.  Engng..  in  press. 

(6)  "On  the  Variational  Structure  of  Assumed  Strain  Methods,”  with  T.J.R.  Hughes.  J.  Appl. 
Mechanics,  in  press. 

(7)  "An  Analysis  of  a  New  Class  of  Integration  Algorithms  for  Elasioplastic  Constitutive  Rela¬ 
tions,"  with  M.  Ortiz,  Int.  J.  Num.  Meth.  Engn..  in  press. 

(8)  "A  3-Dimenslonai  Finite  Deformation  Viscoelastic  Model  Accounting  for  Damage  ElTects.” 
with  R.L.  Taylor,  Report  No.  UCB/SESM-85/02.  Department  of  Civil  Engineering.  Univer¬ 
sity  of  California,  Berkeley.  1985. 

(9)  "Bending  and  Membrane  Elements  for  the  Analysis  of  Thick  and  Thin  Shells, 

(10)  "Augmented  Lagrangian  Formulations  for  the  Finite  Element  Solution  of  Contact  Prob¬ 
lems."  with  P.  Wriggers  and  R.L.  Taylor,  NUMETA  Conference,  Swansea,  Wales,  January 
1985 

(11)  "Rate  Constitutive  Equations;  Are  they  Necessary?",  with  M.  Ortiz,  K.S.  Pister  and  R.L. 
Taylor.  ASME  Winter  Annual  Meeting.  New  Orleans  1984. 


